Neutrinos and the Highest Energy Cosmic Rays 



O 
O 

o 

(N 



> 

o 
o 
o 

o 



% 



A. Letessier-Selvon^* 

'^Laboratoire de Physique Nucleaire et des Hautes Energies, 

IN2P3-CNRS, Universites Paris 6 & 7, 

4 place Jussieu, Tour 33 RdC, 75252 Paris Cedex 05 France 

Observation of Ultra High Energy Cosmic Rays (UHECR) -whose energy exceeds 10^*^ eV- is still a puzzle for 
modern astrophysics. The transfer of more than 16 Joules to a microscopic particle can hardly be achieved, even 
in the most powerful cosmic accelerators such as AGN's, GRB's or FR-II radio galaxy lobes. Potential sources 
must also lie within 100 Mpc of the Earth as the interaction length of protons, nuclei or photons is less than 
lOMpc. However no visible counterpart of those sources has been observed. Calling upon new physics such as 
Topological Defect interactions or Super Massive Relic Particle decays is therefore very tempting, but such objects 
are yet to be proven to exist. Due to the very low flux of UHECR only very large dedicated experiments, such as 
the Auger observatories, will allow to shed some light on the origin of those cosmic rays. In this quest neutrinos, 
if they can be detected, are an invaluable messengers of the nature of the sources. 



1. INTRODUCTION 

The cosmic ray spectrum is now proved [|l|,|[ to 
extend beyond 10^°eV. To be observed on Earth 
with such energies, particles must be produced or 
accelerated in the Universe with energy near or 
above lO^^eV. Conventional acceleration mecha- 
nisms in astrophysical objects can only reach this 
requirement by stretching to the limit their avail- 
able parameter space, making such scenarios un- 
likely to explain the origin of UHECR. Alterna- 
tive hypothese involving collapse of Topological 
Defects (TD) or decay of Super Massive Relic 
Particles (SMRP) are well suited to produce par- 
ticles above 10^°eV but need a proof of existence. 

Transport, from the source to Earth, is also 
an issue. At those extreme energies the Cos- 
mic Microwave Background Radiation makes the 
Universe essentially opaque to protons, nuclei 
and photons which suffer energy losses from pion 
photo production, photo-disintegration or pair 
production. These processes led Greisen, Zat- 
sepin and Kuzmin[g| to predict a spectral cut- 
off in the cosmic ray spectrum around 5 x lO^^eV, 
the GZK cutoff. The available data although still 
very scarce do not support the existence of such 
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a cutoff. Therefore the sources are either close by 
and locally over dense for the cutoff not to show, 
or new physics prevents the UHECR from the ex- 
pected energy losses against the CMB photons. 

The following will briefly develop the argu- 
ments mentioned in this introduction. Interested 
readers should consult the numerous reviews de- 
voted to the subject H-0] for more details. 



2. OBSERVATIONS 

The differential spectrum of cosmic ray flux|7] 
as a function of energy is shown on Figure 1. 
Integrated fluxes above three energy values are 
indicated: 1 particle/m^-second above 1 TeV, 1 
particle/m^-year above 10 PeV, 1 particle/km^- 
year above 10 EeV. The energy spectrum is sur- 
prisingly regular in shape. From the GeV energies 
to the GZK cutoff, it can be represented simply 
by three power-law segments interrupted by two 
breaks, the so-called "knee" and "ankle" . 

Figure g is a zoom on the highest energy part of 
the total spectrum where only the latest AGASA 
data|^ is displayed. On this figure, the energy 
spectrum is multiplied by E^ so that the part be- 
low the EeV energies becomes flat. Comparing 
the data points and the dashed line one has a 



X 10 



Fluxes of Cosmic Roys 



(1 particle per m^-second) 






Knee 

(1 particle per m^-year) 



Ankle 

(1 particle per km^-year) 




Energy (eV) 



Figure 1. The cosmic rays spectrum [Q 



clear view of what can be expected from a cos- 
mological (uniform) distribution of conventional 
sources and what is observed. 

The cutoff, that would be expected if the 
sources were cosmologically and uniformly dis- 
tributed and if the observed cosmic rays had no 
exotic propagation or interaction properties, is 
not present in the observed data. 

In the search for potential sources, one looks for 
correlations of the UHECR arrival directions with 
the distribution of matter within a few tens of 
Mpc. Such an analysis was done by the AGASA 
experiment for the highest energy range |^. No 
convincing deviation from isotropy was found. 

If the sources of UHECR are nearby astrophys- 
ical objects and if, as expected, they are in small 
numbers, a selection of the events with the largest 
magnetic rigidity would combine into multiplets 
(cluster of events whose error boxes overlap). 

Figure || shows the subsample of events in 
the AGASA catalog with energies in excess of 
100 EeV (squares) and in the range 40-100 EeV 
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Figure 2. Highest energy region of the cosmic ray 
spectrum as observed by the AGASAB detector. 



(circles). One can see that there are three dou- 
blets and one triplet. The chance probability of 
having as many multiplets as observed with a 
uniform distribution is estimated to be less than 

1% i- 

The non uniform sky coverage -all present de- 
tectors are in the northern hemisphere- and the 
small statistics available make anisotropy studies 
difficult. The Auger observatories are designed 
with full sky coverage and large detection areas 
to overcome these difficulties. 

3. TRANSPORT AND PRODUCTION 

Today's understanding of the phenomena re- 
sponsible for the production of UHECR is still 
limited. One distinguishes two classes of pro- 
cesses: the "Top-Down" and "Bottom-Up" sce- 
narios. In Top-Down scenarios, the cosmic ray 
is a decay products of a super-massive particle. 
Such particles with masses exceeding lO^^eV are 
either meta-stable relics of some primordial field 
or a GUT gauge boson produced by the radia- 
tion or collapse of topological defects. In the 
Bottom-Up scenarios, the energy is transferred to 
a charged particle at rest through its electromag- 
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Figure 3. Arrival directions (galactic coordinates) 
of cosmic rays with E >40 EeV, AGASA.H 



netic interactions. This classical approach does 
not require new physics. 

At energies above 10 EeV and except for neutri- 
nos, the Universe is not transparent to ordinary 
stable particles on scales larger than about 100 
Mpc. Regardless of their nature, cosmic rays lose 
energy in their interaction with the various pho- 
ton backgrounds, dominantly the copious Cosmic 
Microwave Background (CMB) but also Infra- 
Red and Radio. The absence of prominent visible 
astrophysical objects in the direction of the ob- 
served highest energy cosmic rays together with 
this distance limitation adds severe constraints on 
the "classical" Bottom-Up picture. 

3.1. GZK cutoff 

The energy at which the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff takes place is given by the 
threshold for pion photo-production in the pro- 
ton CMB-photon collisions. For an average CMB 
photon (IQ-^gV), one obtains Eth = TxlO^^eV. 
The interaction length can be estimated from the 
pion photo-production cross section and the CMB 
density : L — {crp)~^ ~ 6 Mpc. 

A recent Monte Carlo [gj calculation, includ- 
ing red shift, pair production and pion photo- 
production losses, is shown on Figure 0. The loss 
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length xioas is defined as xios 
100 EeV photo-production processes are domi- 
nant and the loss length falls below 13 Mpc. 

For nuclei, the situation is usually worse. They 
photo-disintegrate in the CMB and infrared radi- 
ations losing on average 3 to 4 nucleons per Mpc 
when their energy exceeds 2 x lO^^eV to 2 x lO^'^eV 
depending on the IR background density value. 
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Figure 4. Top : Loss length of a proton from M 
Bottom : ratios with other calculations. 



Top-Down production mechanisms predict 
that, at the source, photons and neutrinos domi- 
nate over ordinary hadrons by about a factor four 
to tenPpd|]. An observed dominance of gammas 
or neutrinos in the supra-GZK range would then 
be an inescapable signature of a super-heavy par- 
ticle decay or TD interaction. High energy pho- 
tons traveling through the Universe produce e^e~ 
pairs when colliding with the Infra-Red/Optical 
(IR/0), CMB or Universal Radio Background 
(URB) photons. As can be seen on Figure pi the 
attenuation length gets below 10 Mpc for photon 
energies between SxlO^'^eV and lO^'^eV. In this 
energy range the Universe is opaque to photons 
on cosmological scales. 

Once the photon has converted, the e'^e^ pair 
will in turn produce photons mostly via Inverse 
Compton Scattering (ICS). Those two dominant 
processes are responsible for the production of 
electromagnetic (EM) cascades. On Figure ^ one 
sees that, at the pair production threshold on 
the CMB photons (2xlO^''eV), conversion occurs 
on distances of about 10 kpc (a thousand times 
smaller than for protons at GZK energies) while 
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Figure 5. Attenuation length of photons, protons 
and iron. Double pair production (not shown) 
limits the photon attenuation length to about 
lOOMpc above lO^^eV.pOl 



subsequent ICS of electrons on the CMB in the 
Thomson regime will occur on even smaller scales 
(1 kpc). 

As a consequence, most photons of ultra high 
energy will produce, through successive collisions 
on the various photon backgrounds (URB, CMB, 
IR/0), lower and lower energy cascades and pile 
up in the form of a diffuse photon background 
below 10^^ eV with a typical power law spectrum 
of index a = 1.5. This is a very important fact 
as measurements of the diffuse gamma ray back- 
ground in the 10^-lO^^eV range done for example 
by EGRET ||l^] will impose limits on the photon 
production fluxes of Top-Down mechanisms and 
consequently on the abundance of topological de- 
fects or relic super-heavy particles. 

Neutrinos are the only known particles that 
can travel through space unaffected even on large 
distances, carrying intact the properties of the 
source to the observer. They may prove to be an 
unambiguous signature of the new physics under- 
lying the production mechanisms. 



3.2. Bottom-Up acceleration 

In the conventional acceleration scenarios one 
distinguishes two types of mechanisms : 

• Direct acceleration by very high electric 
fields in or near very compact objects. This 
does not naturally provide a power-law 
spectrum. 

• Diffusive shock acceleration in all systems 
where shock waves are present. This sta- 
tistical acceleration, known as the Fermi 
mechanism, naturally provides a power-law 
spectrum. 

Hillas has shown ||l^ that irrespective of the de- 
tails of the acceleration mechanisms, the maxi- 
mum energy of a particle of charge Ze within a 
given site of size R is: 
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where B is the magnetic field inside the accel- 
eration volume and /3 the velocity of the shock 
wave or the efficiency of the acceleration mech- 
anism. This condition is nicely represented by 
the Hillas diagram shown in Figure 0. Inspecting 
this diagram one sees that only a few astrophys- 
ical sources satisfy the necessary condition given 
by Eq. ffl). Let us briefly review them : 

Pulsars : From a dimensional analysis, one 
obtains up to lO'^^ volts for the potential drop 
in a rotating magnetic pulsar. However the high 
radiation density in the vicinity of the pulsar will 
produce e^^e^ pairs which reduce the potential 
drop down to values of about 10^'^eV. A different 
mechanisms involving Fe nuclei acceleration by 
relativistic MHD winds has been put forward |1J] . 
But details of the effectiveness of this mechanism 
still need to be demonstrated. 

AGN cores and jets : Blast waves in AGN 
jets could in principle lead to a maximum en- 
ergy of a few tens of EeV jl^] and similarly for 
AGN cores. However those maxima are unlikely 
to be achieved under realistic conditions due to 
the interaction of the accelerated protons with 
very high radiation fields in and around the cen- 
tral engine of an AGN. 
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Figure 6. Size and magnetic field strength of pos- 
sible acceleration sites. Objects below the diago- 
nal lines cannot accelerate the corresponding ele- 
ments above 10^°eV or lO^^eV. 



However in such a framework the UHECR spec- 
trum should clearly show the GZK cut-off. 

3.3. Top-Down production 

One way to overcome problems related to the 
acceleration of UHECR and the invisibility of 
their sources is to introduce a new unstable or 
meta-stable super-massive X-particle. The decay 
of this X-particle produces, among other things, 
quarks and leptons, resulting in a large cascade 
of energetic photons, neutrinos and light leptons 
with a small fraction of protons and neutrons, 
part of which becomes the UHECR. For this sce- 
nario to be observable three conditions must be 
met: 

• The decay must have occurred recently 
since the decay products must have trav- 
eled less than about 100 Mpc because of 
the attenuation processes discussed above. 

• The mass of this new particle must be well 
above the observed highest energy (100 EeV 
range), a hypothesis well satisfied by Grand 
Unification Theories (GUT) whose scale is 
around lO^^-lO^^eV. 



FR-II radio galjLxies : Radio-loud quasars 
are characterized by a very powerful central en- 
gine ejecting matter along thin extended jets. At 
the ends of those jets, the so-called hot spots, 
the relativistic shock wave is believed to be able 
to accelerate particles up to ZeV energies. FR- 
n galaxies seem the best potential astrophysical 
source of UHECR ||I^ . Unfortunately, no nearby 
(less than 100 Mpc) object of this type is visible 
in the direction of the observed highest energy 
events. 

Gamma Ray Burst : Gamma ray bursters 
(GRB) are intense sources of gamma rays. The 
most favored GRB emission model is the "ex- 
panding fireball model" where one assumes that 
a large fireball, as it expands, becomes optically 
thin hence emitting a sudden burst of gamma 
rays. The observation of afterglows allowed mea- 
surement of the red shift of the GRBs hence con- 
firming their cosmological origin. GRB can be 
shown to accelerate protons up to 10^"eV|17|. 



• The ratio of the volume density of this par- 
ticle to its decay time must be compatible 
with the observed flux of UHECR. 

The X-particles may be produced by way of two 
distinct mechanisms: 

• Radiation, interaction or collapse of Topo- 
logical Defects (TD), producing X-particles 
that decay instantly. In those models the 
TD are leftovers from the GUT symmetry- 
breaking phase transition in the very early 
universe. Quantitative predictions of the 
TD density that survives a possible infla- 
tionary phase rely on a large number of the- 
oretical hypotheses. Therefore they cannot 
be taken at face value, although the ex- 
perimental observation of large differences 
could certainly be interpreted as the signa- 
ture of new effects. 

• Super-massive meta-stable relic particles 
from some primordial quantum field, pro- 



duced after the now commonly accepted in- 
flationary stage of our Universe. However 
the ratio of their hfetime to the age of the 
universe requires a fine tuning with their 
relative abundance. It is worth noting that 
in some of those scenarios the relic particles 
may also act as non-thermal Dark Matter. 

In all conceivable Top-Down scenarios, photons 
and neutrinos dominate at the end of the hadronic 
cascade. This is the important distinction from 
the conventional acceleration mechanisms. 
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Figure 7. Horizontal shower development 



4. THE AUGER DETECTOR 

Large area ground based detectors do not ob- 
serve the incident cosmic ray directly but the Ex- 
tensive Air Shower (EAS), a very large cascade 
of particles, they generate in the atmosphere. All 
experiments aim to measure, as accurately as pos- 
sible, the direction of the primary cosmic ray, its 
energy and its nature. There are two major tech- 
niques used. One is to build a ground array of sen- 
sors spread over a large area, to sample the EAS 
particle densities on the ground. The other con- 
sists in studying the longitudinal development of 
the EAS by detecting the fluorescence light emit- 
ted by the Nitrogen molecules which are excited 
by the EAS secondaries. 

The Auger Observatoriesn combine both tech- 
niques, with construction starting in the fall of 
2000. Once completed in 2006, they will be cov- 
ering two sites, one in the southern hemisphere 
(Argentina) and one in the north (Utah, USA). 
The surface of each site, 3000 km'^, will provide 
statistics of a few tens of events per year above 
100 EeV. The detector is designed to be fully effl- 
cient for showers above 10 EeV, with a duty-cycle 
of 100%. Each station of the ground array is a 
cylindrical Cerenkov tank of 10 m^ surface and 
1.2 m height flUed with filtered water. Because 
of the size of the array, the stations have to work 
in a stand-alone mode: they are powered by solar 
panels and batteries, communication is wireless 
and timing is provided by the GPS satellites. 

The fluorescence telescopes use photo-tubes 
with a field of view of 1.5°. Each telescope 
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sees an angle of about 30 x 30 degrees. On the 
southern site, three eyes (7 telescopes each) will 
be installed at the periphery of the array and 
one (12 telescopes) in the middle, in order for 
the whole array to be visible by at least one of 
the telescopes. In the hybrid mode (10% of the 
events) , the detector is expected to have on aver- 
age 10% energy resolution and an angular preci- 
sion of about 0.3°. 

With a total aperture of 14000 km^sr (both 
sites), the Auger Observatory should detect every 
year of the order of 10000 events above 10 EeV 
and 100 above 100 EeV. 

5. NEUTRINOS 

Although both neutrinos and photons domi- 
nate the particle fluxes in Top-Down models, only 
neutrinos are perfect probes of the characteristics 
of the sources. High energy neutrinos can also be 
produced in Bottom Up scenarios as secondaries 
of hadronic interaction. If AGNs, Radio Galaxy 
lobes or GRBs are UHECR sources they should 
produce a substantial flux of neutrinos. Finally, 
along their path in the universe, hadrons will also 
produce neutrinos via the pion photo -production 
processes, the GZK neutrinos. 

The first study on the detection of UHE neutri- 
nos with the Auger detector were done by [Q^ . 
The UHE neutrinos may be detected and distin- 
guished from ordinary hadrons by the shape of 
the horizontal EAS they produce. At large angles, 
above 60°, hadronic showers have their electro- 
magnetic part extinguished and only high energy 
niuons survive. Therefore the shower front is very 
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Figure 8. Neutrino fluxes from various 

sources [Q, dotted lines speculative, dashed 
probable, solid certain. The 2 top thick solid 
lines represent the Auger sensitivity (0.3 event 
per year) p^]. 



flat (radius of curvature is larger than 100 km), 
and very narrow (less than 50 ns). Neutrinos in- 
teracting deeply in the atmosphere will start a 
shower above the detector which will appear as 
a "normal" shower, with a curved front (a few 
km), a large electromagnetic component, and a 
wider signal (a few microseconds) [see Figure g. 
With such important differences and if the fluxes 
are high enough, neutrinos will be identified and 
detected . 

Figure shows the expected fluxes from a 
model calculation by Protheroe p8|. The sen- 
sitivity limit of the Auger detector dcflncd as 0.3 
events per year is also shown. Although each site 
of the Auger observatory reaches 10 km^ water 
equivalent of target mass, only the models classi- 
fled as speculative by the author are expected to 
yield a detectable signal. 

6. CONCLUSIONS 

The composition, the shape of the energy spec- 
trum, and the distribution of arrivals of UHECR 
will prove to be powerful tools to distinguish be- 
tween the different production scenarios. 

If UHECR are hadrons accelerated by Bottom- 
Up mechanisms, they should point back to their 
sources, with visible counterparts. 

For Top-Down mechanisms and above 100 EeV, 
one should observe a flux of photons and neutri- 
nos as the photon absorption length increases (up 
to a few 100 Mpc). Below 100 EeV the spectrum 
shape and composition will depend on the charac- 
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teristic distance between TD interactions or relic 
particle decays and Earth, the proton attenuation 
length and the photon absorption length. 
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